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We study the specific heat scaling function of superfluids confined in cubic geometry and in 
parallel-plate (film) geometry with open boundary conditions (BC) along the finite dimensions 
using Monte Carlo simulation. For the case of cubic geometry for the superfluid order parameter 
we apply open BC in all three directions. We also calculated the specific heat scaling function for 
the parallel-plate confinement using open BC along the finite dimension and periodic BC along 
the other two dimensions and we find it to be very close to the earlier calculated using Dirichlet 
instead of open BC. We find that the specific heat scaling function is significantly different for 
the two different geometries studied. In addition, we generally find that the scaling function for a 
fixed given geometry when calculated with open BC is quite close to that calculated with Dirichlet 
BC, while when calculated with periodic BC is quite different. Our results for both the scaling 
functions obtained for the parallel-plate geometry and for cubic geometry with open BC along 
the finite dimensions are in very good agreement with the recent very high quality experimental 
measurements with no free parameters. 



I. INTRODUCTION 

Thermodynamic quantities such as the specific heat 
become non-analytic at a critical point associated with a 
second order phase transition. For a finite (or confined) 
system with a finite dimension such as a film character- 
ized by a length L, close enough to the critical point such 
that the correlation length becomes comparable or larger 
than L, such thermodynamic quantities are significantly 
altered; the reason is that the degrees of freedom of the 
system arc correlated to each other over the entire sys- 
tem. Examples of such confined systems are (a) a film 
of thickness L where the system is confined in one spa- 
tial dimension, (b) a bar-like geometry with cross-section 
Lx L and infinite length (such a pore or a wire) where the 
system is confined in its two spatial dimensions or (c) a 
cubic geometry of size I? where the system is finite in all 
three dimensions. For any thermodymamic observable 
we can define a system specific dimensionless quantity, a 
scaling functionjj]; for example, in the case of the specific 
heat near the critical point and for sufficiently large L we 
may define the following scaling function: 



f(x = tL 1 ^) 



c(t,L) - c(0,oo) 
c(0, L) — c(0, oo) : 



(1) 



where t = T/T\ — l is the reduced temperature and c(t, L) 
is the specific heat for the case of the system confined 
within a finite length L. For a given case of confining 
geometry and given the condition which the order pa- 
rameter satisfies at the boundaries of the system, as L 
approaches infinity and t — > 0, the scaling function de- 
pends only on the value of the combination x = tL~ , 
namely on the length L in units of the correlation length 
£(t) ~ t~". A dimensionless function such as f(x) de- 
fined by Eq. ^ can be thought of as a universal scaling 



function for the specific heat for a well-defined confining 
geometry. In other words, the scaling function does not 
depend on the microscopic details, but only depends on 
the nature of the universality class of the system, the 
confining geometry and the boundary conditions which 
are felt by the order parameter. 

In this limit (t — ► and L — > oo) the scaling function 
is different for the three different cases mentioned previ- 
ously for the following reason: For a fixed value of x << I 
and for any large value of L there is always a sufficiently 
small value of t satisfying the condition where the correla- 
tion length is much larger than L. However, in this limit 
the case (a) is the case of a 2D plane, the case (b) is the 
case of a ID line and the case (c) is the zero-dimensional 
case. Thus the value f(x) for sufficiently small values of 
|x| should be very different for these three geometries. 

Though earlier experiments on superfluid helium films 
of finite thickness seemed to confirm the validity of the 
finite-size-scaling(FSS), there were later experiments U 
U where it was shown that the superfluid density of thick 
helium films does not satisfy FSS when the expected val- 
ues of critical exponents were used. Similarly, in mea- 
surements of the specific heat of helium in finite geome- 
tries, other than the expected values for the critical ex- 
ponents were found 

More recent experiments in microgravity 
environment Q as well as earth bound experiments^, H] 
are consistent with scaling and they have determined 
the specific heat scaling function for the parallel plate 
(film) geometry (case (a)) and they are in reasonable 
agreement with the scaling function as was predicted 
by Monte Carlo simulations 0, ^} and renormalization 
group techniques 0. While the specific heat scaling 
function for case (b) confinement has been theoretically 
determined|l2j and it was found to be significantly 



2 



suppressed relative to the case (a) there are so-far 
no experimental data to compare. More recently, the 
specific heat scaling function for the case (c) has been 
experimentally determined p^. Il^j . 

The main goal of this paper is to present the results 
of our Monte Carlo simulations to determine the spe- 
cific heat scaling function for cubes with open boundary 
conditions (BC) in all three directions (confining case 
(c)). In this case the scaling function characterizes the 
zero-dimensional to three-dimensional transition. Our 
results for the scaling function are compared to the very 
recently obtained experimental results for specific heat 
scaling function in the case of cubic confinement p^. Il4| . 
We find satisfactory agreement with no free parameters. 
In addition, we present results for the specific heat scal- 
ing function for the parallel plate geometry on lattices of 
size L\ x L2 x L with L± = L2 » L where we have ap- 
plied periodic BC along the Li^-directions and open BC 
along the film-thickness direction of size L. The latter 
case was carried out in order to compare the results for 
Dirichlet BC (on the top and bottom of the film) obtained 
earlier 0| . In Refs it was found that while the 

results with periodic BC along the film-thickness direc- 
tion were very different from those obtained with Dirich- 
let BC, the results obtained with Dirichlet BC fit the 
experimental results with no free parameter. In this pa- 
per we find that the scaling function obtained with open 
BC along the finite dimension is close to that obtained 
with Dirichlet and also fits reasonably well the experi- 
mental results obtained by the so-called Confined Helium 
Experiment^ (CHEX). 



II. MONTE CARLO CALCULATION 

We have performed a numerical study of the scaling 
behavior of the specific heat of 4 He in a cubic and in a 
film geometry at temperatures close to the critical tem- 
perature T\. The superfluid transition of liquid 4 He be- 
longs to the universality class of the three-dimensional 
x — y model, thus, we are going to use this model to com- 
pute the specific heat at temperatures near T\ using the 
cluster Monte-Carlo method |T|j . The x — y model on a 
lattice is defined as 



expression 



H = 



(i,j) 



J ' 



(2) 



where the summation is over all nearest neighbors, s = 
(cos 9, sin 9) is a two-component vector which is con- 
strained to be on the unit circle and J sets the energy 
scale. 

We define the energy density of our model as follows: 



E={e)=Z-UY d S i -s j 

\(i,j) 



(3) 



where V = L 3 for the cubes and V = HL 2 for the film 
geometry. We have calculated the specific heat using the 



c = VT-\{e 2 )-{ef) 



(4) 



The above thermal averages denoted by the angular 
brackets are computed according to 



(O) = Z~ x / JJ ddi 0[9] exp(-pH), 



(5) 



where H = Hj J the energy in units of J and (3 = J/T. 
O[0] denotes the dependence of the physical observable 
O on the configuration and the partition function 

Z is given by 



Z = J\\d9 % exp(-pn). 



(6) 



We computed the specific heat c(T, L) of the x — y 
model as a function of temperature on several cubic lat- 
tices L 3 (with L = 20,30,40,50 ). Open (free) bound- 
ary conditions were applied in all directions, namely the 
spins on the surface of the cube are free to take any value. 
These spins interact only with the 5 nearest neighbors, 
one in the interior and 4 on the surface of the cube and 
there is one missing neighbor. We have also calculated 
the specific heat scaling function fi(x) (to be defined in 
the following section) for the case of the parallel plate ge- 
ometry L\ x Li x L (£1,2 >> L) using periodic boundary 
conditions along the long directions of the film and open 
BC along the thickness direction L. For this case we need 
to take the limit L12 — > 00 first; in Ref. |10) it was found 
that using L\ = L2 = 5L was large enough, in the sense 
that systematic errors due to the finite-size effects from 
the fact that L± t 2 are not infinite are smaller than the 
statistical errors for realistic computational time scales. 
The present simulations for films were done on lattices 
60 x 60 x 12, 70 x 70 x 14, and 80 x 80 x 16. 

We used the Monte Carlo method and in particular 
Wolff's cluster algorithm[l5| . We carried out of the order 
of 30,000 MC steps for thcrmalization to obtain equilib- 
rium configurations. We made of the order of 10,000- 
50,000 measurements allowing 500 MC steps between 
successive measurements to obtain statistically uncorre- 
cted configurations. 



III. SCALING FUNCTIONS 

The main goal of this paper is to present a calculation 
of the specific heat scaling function for the case of cubic 
confinement and open BC. In this calculation we have 
used open BC in all three directions of the cube. We 
found that using open BC the results for the specific heat 
scaling were very close to those obtained with Dirichlet 
BC. This will be demonstrated in the next section where 
we compare the previously published results 0, llOj for 
films with Dirichlet BC and results reported in this paper 
for films with open BC. 
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FIG. 1: The scaling function —f(x) denned by Eq.Q 



One can imagine a number of different scaling func- 
tions for the specific heat. Any dimensionless combina- 
tion such as the ratio given by Eq.^can be used as a scal- 
ing function. However, the various experimental groups 
have extracted two scaling functions, the so-called fi(x) 
and f 2 {x) with x = tL~ . These scaling functions are 
defined as follows: 

c{t,L)-c{t Q ,o Q ) = L<*l v f 1 {tl}/ v ) (7) 

[c(t,oo)-c(t,L)]t a = f 2 (tL 1 ^) (8) 

We limit our goal to calculate the specific heat scaling 
function for the confined geometry and not the bulk crit- 
ical exponents nor critical amplitude ratios. We take the 
values for the bulk critical exponents and the universal 
amplitude ratios as determined experimentally^^. Pre- 
vious MC work such as the work of Ref. shows that 
the critical exponents are within error bars from the ex- 
perimental values. Our approach to use the experimen- 
tally determined values of the critical exponents and am- 
plitude ratios and to determine the scaling function by 
applying FSS on the calculated c(t, L), has no fitting pa- 
rameters and this allows no ambiguity. Therefore we use 
v = 0.6709 as obtained from accurate experiments^^ 
such as the so-called Lambda Point Experiment (LPE), 
an experiment in microgravity environment. The hyper- 
scaling relation a = 2 — Sv yields a/v = —0.0189, and 
the correlation length = becomes equal to 

the system size L at the reduced temperature to, i.e., 
to = {iv/L) 1 / u with Co = 0.498. 

In order to find the universal function f(x) defined by 
Eq. ^ we need to know c(0, oo). We use the bulk values 
for c(0, oo) obtained by studying the finite-size scaling of 
the specific heat of cubes with periodic BC^j. I n Fig- 01 
the scaling function —f(x) obtained for cubes with open 
BC in all three directions is compared to that obtained 
with periodic BC|l7|. 



The scaling function fi (x) (Eq. 0) can be calculated 
using our calculated c(t, L) and 

c(io,co) = c(0 ) oo) + c^ a , (9) 

where the values of c(0, oo) and for the x — y model 
are obtained from reference|l7|. 

In order to calculate the universal function fo(x) 
(Eq. , we need to know the bulk specific heat c(t, oo) 
also. Since we are restricting ourselves to the critical 
region we may write the following 

c(t >0,oo) = c(0,oo) + c+r a , (10) 
c(i <0,oo) = c(0,oo)+ct/r\t\- a , (11) 

r = 4 (12) 
c i 

where r is the universal amplitude ratio and it is most 
accurately determined exp er iment all v 1 1 8l [l9j from the 
critical properties of bulk helium to be r = 1.053(2)^8|. 
Inserting Eqs. E3 and EU into Eq. GO we obtain 

f 2 {tL L / v ) = [c(0,oo)-c(t,L)]t a + c+, t>0, (13) 
h{tL l/v ) = [c(0,oo)-c(t,L)]\t\ a + c+/r t < 0. (14) 

which can be calculated by using our computed c(t, L) 
and the values of c(0, oo) and cf from Ref. [l7| . 

IV. FILMS WITH OPEN BOUNDARY 
CONDITIONS 

In Ref. 0, the specific heat scaling function for a 
parallel plate geometry on lattices of size L\ x L 2 x L with 
L\ = L 2 >> L was calculated. In Ref. 9, 10] periodic BC 
along the Li.2-directions and staggered (Dirichlet) BC or 
periodic along the film-thickness direction of size L were 
applied. It was found that while the calculated scaling 
function for the parallel-plate geometry using periodic 
BC along all three directions was very different from that 
obtained with Dirichlet BC along the top and bottom of 
the plate and periodic BC along the other two long di- 
rections, the latter scaling function fits the experimental 
results with no free parameter. This was explained on 
the basis that physically the order parameter along the 
film thickness vanishes at the boundaries of the film and 
therefore Dirichlet BC are more appropriate. 

In this paper we have used open BC along the top and 
the bottom of the plate, instead of Dirichlet, and periodic 
BC along the two long directions of the plate. Since the 
film terminates on the top and on the bottom surface, 
for the pseudospins which belong to these two surfaces 
(in language of the x — y model) there is no neighboring 
spins beyond the top and the bottom surface plane of 
the plate. Therefore, even if we use open (free) boundary 
conditions this termination acts as "a zero order param- 
eter constraint" beyond the top and the bottom of the 
plate. This implies that these two BC, namely staggered 
BC and open BC, are very similar for thick enough films. 
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In order to make a direct comparison of our calculated 
fi(x) to the experimental f%(x), we express all lattice 
units in physical units using the following equation^ 

fl{x)\phys = A/i(x) {lattice, (15) 

A =V^( a /A)-Z (16) 

where V m is the molar volume of liquid helium at the 
lambda point and saturated vapor pressure, fcs is Boltz- 
mann's constant and a the lattice spacing in the x — y 
model required to make contact with the critical behav- 
ior of the correlation length in helium. This prefactor 
A = 15.02 J I (K mol) and it was determined in Ref.Q- 

In Fig. |2 we compare the results for j\ for the case of 
films obtained with open BC along the direction of the 
film thickness to those obtained earlierjij El and to the 
experimental results 0, 0. It is clear that within error 
bars our results for the specific heat scaling function are 
the same for both cases of BC. 
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FIG. 2: Film geometry: The computed universal function 
fi(x) with open BC (solid circles) is compared to the previ- 
ously calculated scaling function using Dirichlet BC0(data 
shown as stars) and periodic BC (shown as plus signs) and 
the experimental results of Lipa el «Z.|g (open circles) and 
those of Mehta el aZ.|g (open triangles). 

While the scaling function is sensitive to boundary con- 
ditions, this indicates that it is hard to distinguish Dirich- 
let from open BC for the specific heat scaling function. 
We feel that when we use physical BC the agreement 
between the theoretical results for the specific heat scal- 
ing function and the experimental results is quite good 
taking into consideration the fact that there is no free 
parameter. 



V. CUBIC CONFINEMENT 

In this section we present the results for the scaling 
functions fi(x) and fiix) obtained for cubes of size L 3 
with L = 20,30,40,50 using open and periodic BC in 



all three directions. As was shown in the previous sec- 
tion open (free) BC are similar to using Dirichlet BC and 
they both express the physical condition imposed by the 
confinement or the termination of the system. In Fig. [3] 
we compare the scaling function fi (x) obtained for open 
BC with that obtained for periodic BC^%. Notice the 
suppression of /i(x) when calculated with open BC rel- 
ative to the case of periodic BC. This is similar to the 
case of the parallel-plate geometry (Fig.0). The scaling 
functions /i(x), however, are very different for cubic and 
parallel-plate geometry. Notice, for instance, that for the 
case of cubic confinement with open BC fi(x) is negative 
for all values of x something very different of what hap- 
pens for any of the calculated of the experimental scaling 
functions for parallel-plate confinement. 

1.5 | . , . , . , 1 
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FIG. 3: The scaling function fi(x) obtained for cubes of size 
L with open (solid circles) and that obtained for cubes with 
periodic BC (open circles) are compared. 

In Fig. 0] we give the results of our present Monte Carlo 
calculation of the function fi (x) for cubes with open BC 
using Eqs 1131141 Fig. [3] shows the results of our calcula- 
tion with periodic boundary conditions. Fig. [S] compares 
the scaling function jiix) obtained for open BC and for 
periodic BC. Notice the qualitatively different behavior 
for the same scaling function for the same geometry but 
different boundary conditions. 

Experimentally the universal scaling function f-iix) for 
cubic confinement has just become available [FJ |l4|. In 
order to make a direct comparison of our calculated scal- 
ing function f%(x) to the experimentally determined, we 
express all lattice units in physical units. The prefactor 
is the same as in the case of the function fx (x) : 

h{ x )\phys = Xf2(x)\lattice, (17) 

where A is the constant given in the previous section by 
Eg. 1161 and its numerical value is A — 15. 02 J/(K mol). 

In Fig. El/2 (x) obtained from our MC calculation is 
compared with the experimental data0,0|. The agree- 
ment between the scaling function calculated with open 
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FIG. 4: The universal function f2(x) for cubes of size L 3 with 
open boundary conditions. 
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FIG. 5: The universal function f2(x) for cubes of size L with 
periodic boundary conditions 
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FIG. 6: The computed universal function ji(x) for open and 
periodic BC and for cubes is compared to the experimental 
restdtsplltH. 



BC and experiment is quite satisfactory considering the 
fact that there is no free parameters. 



VI. CONCLUSIONS 

In this paper we have used the x — y model which de- 
scribes the fluctuations of the superfluid order parameter 
near the critical point to calculate the scaling functions 
associated with the specific heat for the case where the 
superfluid is confined in a cubic geometry and in parallel- 
plate geometry. Both in the theoretical calculations and 
in the experiments, the region very near the superfluid 
transition is probed such that the correlation length asso- 
ciated with the superfluid order parameter is of the size 
of the confining length. 

First, we calculated the specific heat scaling function 
for the case of parallel plate confining geometry using 
open boundary conditions along the top and the bottom 
surfaces of the film. Our results are very close to those 
obtained^ with Dirichlet (staggered BC) along the 
top and the bottom surfaces of the film. Both calcula- 
tions are in satisfactory agreement with the experimental 
results [USUI while the results of earlier calculations us- 
ing periodic boundary conditions [jjj were found to dis- 
agree with the experimental scaling function near the su- 
perfluid transition. 

Just recently, experimental measurements 0, \^ on 
superfluid helium confined in cubes became available. 
This prompted us to calculate the heat capacity scaling 
function of superfluids for cubic confinement. When we 
used open boundary conditions in all three directions of 
the cube we find very good agreement between the cal- 
culated and the measured [jjj UJj scaling functions with 
no adjustable parameter. On the contrary, if periodic 
boundary conditions are used at the boundaries of the 
cube, which are unphysical boundary conditions for a 
confined system, there is great disagreement between the 
calculated and the measured specific heat scaling func- 
tions. 
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